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ABSTRACT 


This paver deals witn the refraction of long 
Gravity waves ¿bout a circular island with unifera 
slope below mean water level. The Fermat principle is 
used as a basis to derive an analvtical solution for 
the ray e A refraction analysis is performed to 
determine the refraction factor, phase lag, and wave 
heignt at the island shore. “he relative power distri- 
Durt.on of rays in the far field is determined by con- 
sidering emergent reflected and refracted waves from the 


islande 


n P 
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INTRODUC TATY 


0 ۹۶۳ئ۴ 5۹ تح 7۹5ج جج 237-5 


The Purpose of rt 6 ብ ቤዬዎኞቕዎ፡ "۸6 ٦21 “ene tee ዋጮ. 
on long gravity waves of an island with circular bottom 
MBoOnecOurs and radially Gemstane slope. Lsing the pring 


iples of wave ray theory, a refraction analysis is 


OQ 
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Troe eour Co determino ai T enod fication of 327٦ 


lane waves as they move into shallow water about 6 


ርን 


Toland The refraction anad direct shoaling effects 


E 


poncxneoming long was ry Q Ra e Considered: fhe r = 


7 |ፎር.ሮመር መ refracted rays emerging from tne island 


are considered in order to study the divergent wave 


2 9 7 
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MELOY و و‎ Lar riel p the 1s1and,.,-Inis re 
2--067 236۴۷ تا 6 01000۷۹۰۲۰۰ 730 0ت‎ tO al mOre exacti @2mear 
ment of the subject by Vastano (1966), involving a 


numerical diffraction technique. 
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For centuries man has concerned himself witn the 
effects of islands on surface waves. The early civilize- 
tions of Micronesians and Polynesians utilizec wave pat- 
Berns to mark fishing @rounds end 0 70 
among the islands. Winkler (1901) recorded some of thei 
pechnooguss through interviews Wien native Gees. ren 
(19433 further illustrates this fact and states: 


Every chief and see pilot possessed 
elaborate charts based upon his own exper- 
lence and on knowledge handed down or gain- 
ed from others. These sea charts are made 
with thin strips of the midrib of the leaf- 
lets of the cocomue arranged on a frame 
usuallv rectangular in shape. The knowledge 
they record is indicated oy the arrange- 
ment of the leaf strips relative to one 
another and by the forms given to them 
Dy bending end "“e@ressing a, VCUrVed “Straps in 
dicate the altered direction taken by ocean 
swells when deflected by the presence of 
elves lands Crot Intersections ere, ٣ 
where these meet and tend to produce a confused 
sea. 


More recent interest on the specific subject or wave 
refraction was stimulated by the needs brought about by 
Jorid Wer II amphibious operations. AS a wave approaches 
shallow water, its velocity begins to decrease ap 
ciably at a depth of about one-halr the wave length anc 
the wave height is influenced by the convergence anc div- 
ergence of energy along the shore. As the wave bends 


Lourd eg aand; ህመህር];ሚ 10፡ሀህ“ሠሆ፡... aay De TOMEI 
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the island "wave lee," but often times not to the extent 
that might be expected. This fact was painfully realized 
emring the invasion of sicily when the methods we. DUNS 
and Munk (1947) were used to forecast wave conditions 
MeiOr to tne landing: “Wespite tre use Of these aC e e 
methods, the wave heights an the lee ot the islas s 


greatly underestimated thus severely hampering the landing 


operation. 
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Tae graphical means ot constructing wave ۵ دا‎ 
diagrams was first used by O'Brien (1942). Application 
to the complex ocean bottoms of nature was made by 


E»unson, O'Brien, and Ras ات‎ R OoOo OP and Arthur, Mungo 


H 


saacs (1952). The wave refraction pattern for any bottom 
topography assumes sucn that the waves change direction 
ana are bent in such a manner trat they tend to assume 
the slope of the depth contours. Aerial photographs 
Munk and Arthur, 1952) qualitatively bear this out. 

zt wes learned that the effect of bottom features upon 
EDoVesmecpends noro 626 absolute depen OT وت‎ Duc 
upon the depth relative to the wave length. Grapnical 
eor struecetlon methods by hand are both time consuming and 
tedious. Griswold (1963) has computer programmed the 
graphical construction of wave rays and Momoi (1964) has 
expanded on Griswold's program to include the effects of 


ےہ٣‏ ۹٦۰٠ح‏ ۰ ۲ت ۳" ح3۶2 





in some special cases, analytical solutions to tne 
weve refraction problem have been developed to improve 
s cy which can be erıicieal in,tne’investicatiesse ዚበ: ን 
enercy distribution especially on the wave lee side of is- 
landse Analytical and experimental studies of refraction 
in the region of a caustic have been made by Pierson (1952). 
ከዝ ይ “ከገጋር Of conformal transformation was applied by 
Pocinski (1950). The actual cases treated are restricted 
Dy the matnematical complications that are introduced by 
Beeeercelcular conformal transformation ana by how wel. 
che transform simulates a beach contour. 2 ۳ 
۴۳۰۰٠٠ ٠ط ۶۲۹۹م‎ S principle as e basis Lor an analogy between 
the refraction and the minimal flight path problem. The 
effects of current and depth distribution are considerec. 
Johnson (1947) also discusses the refraction of ocean 
Surface raves by current for deep water waves. 

Munk ana r Artnur (10352) applied ray theory re ea = 
Jate tne wave intensitv alonc a refracted ray. “nis pro- 
cedure has certain advantages over the usual technique 
pcompucing intensity from measurements of the distance 
between rays; nameiy in the areas oft extreme convergence 
or divergence. 

Analytical solutions for tne refraction of waves 
Around an Islanda with ርር ርር? ቴፓርኽርመ ልሙ.” ጋመ Dottom ርር ንከም 


have been obtained by Arthur (1946) through the application 





of Fermat's principle.  Artnur (31951) @eterM=ned Uae cue 
or ont erfeces In the penetration of wave energy مھ‎ 
Ene "wave lee" of an island are generally the result of 


refraction by underwater topography anc 7 nT 


eer eC cy LON « 





፣ ا‎ 





3. General Descrincion of che Island 


ሥ፦ 


toe Island Dotta ቭፕፕ  ፐ ኸኸ፡ሜ፡ሜ፡ኸ፡ኾሥ“ ለሽ ቱ ሥን ኑን ጵ ታንሽንኸሥንሥን፡ጅ፡፪ጅ፡ጅቹጃ፡ጵጅጀዳ፡ጅ.ጂ.፡ gez 
present study is axially symmetrical (i.e., the depth 
contours are concentric circles.) At 5 km. radius fron 
mice sane axis the topeereapny drops vertically fre. 
Eco sea level co a depelweor V0 sAkm., then assunestea rosa 
۲: لا‎ Out to 41 km. redius, Seyond which mé deta 
remains 4.1 km. (see Figure 1). The periods of tke 
LEX cciccted are 2,242.00 200 0178 m»nutes to provide 
ዜር ም 3) in wave lengths within the tsunami Class. 

There was little attempt here to simulate the 
detailed underwater topographical features of any part- 
medular island. However ecomputer program for the 
study of the interaction of monochromatic waves with an 
island of irregular shape is being developed by Reid 
and Vastano (1966). This program permits depth variation 
end the use of a special orthogonal system. 

lt 2S interesting to note that the selected island 
configuration lends itself to a fairly sinple mathemati 
ን ን ን ን 1 7ت ۱7 ۲۲16ح‎ 6 IIE COTE ELGG COC ENVER 
extended slope line intersects the mean sea level at 
.፡›ጸዘታዚፕ[ዚቭዚቭዙጅሞሎሞዚሞቨቪህቨቭህቨህህቨፐ[፡ችሽ፡ችሽችፋሜ፡ችፋችቺቕጅ ت٦‎ 770-٦ 3ج جک‎ mr Ü S ፡፡ rt more robe 
midable problem presents itself with regard to wave 
amplitude determination and is worthy of further in- 


vestigation. 





ዴሬ ASSumpcions 


ln water of constant depth the wave celerity, C, 
for sinusoidal gravity waves oz small amolitude, neglect- 
ing surface tension, is given by the classical relation 
1 በ2, 12232) 


(1) e (g/k) tanh kh , 


where h is the mean depth of water, k = 21[/1 th wave 
number, and g the acceleration due to gravity. A wave 

 ።> (ር 1S shown in i005 mt wnere the wave 6 77 
eee “ortwice tne empl SS በወ rne vertioal distans 
DeLween cwo successive crests, and the period T is the 
time interval between the passage of two consecutive 

e ests- کت‎ ary WavemGlearer folate] teee, in lengen 

the effect of surface tension on wave speed can be ignored 
(Wiegel, 1964). The wave speed is related to wave length 


and period by 
(2) E 


Tnere are two limiting cases of (1) which are 
appropriate to the so called "deebo" anā "shallow" waver 
Lo c DOC CU S or deep water, N/L is very large 
and tanh kh tends toward unity. Under this condition (1) 


and (2) yield 
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pa 
NA 


On the other hand, under the condition h/L= 1/20, 
© comes to a good approximetion 


ee = gn < 


There is, of course, no sh@erp distinction Derween 
"deep" and "shallow" water. The effect of depth of water 
eave characteristics i15 gradual ard waves in any 
finite depth of water are affected by that depth. It is 


interesting to note that weves classified in the shorcer 


eriod Category are not refracted until near shore anrd 


'ህ 


shallow water, while long period waves are refracted 

further offshore in deeper water,  Tsunamis, for example, 

of periods 2 minutes to 1 hour would constantly undergo re- 

E-actlon effects in the open ocean. In a tsunami stuay, 

only waves that travel at a speed of EE are considered. 
As waves propagate into regions of changing deptn, 


y- 


we consider that T remains unchanged while H ean 


(x, 
Es 


Wai Vg 
Caven an initially plene wave of a given se 

wave height moving towards shallow we&ter, the wave height 

Es ceneraloecaàn vary according to the rollowing fectors: 


fects 
ECTS 


th Fh 


(a) Direct shoaling e 
کرک‎ ጋ retracelonis 
ET l ection 





(d) ^" DiIffreoecron 
(e) Energv lost by friccion and vercolegion 


S his study only the G@irst two factors, grege ያለስ 
ever the sloping bathymetry, however, retlection $5 taken 
nio account at the snore. The energy 0 336 7777-7 o 
Peeaerthact.On OCCUrS infareas Of aporeclaple Conversa nici 
Berdivergence of energy ea :d i: shallow water Zornes ue 
as behind breakwaters. Tne cittraction tends to spread 
the wave energy into potential shadow zones. Friction 
Losses on slopes greater than 0.01 are considered nec- 
mictble (Walsh, Reld, ard Bader, 1962). „Percolation 5 
metered to the bottom echeracteristics and 1s considered 
negligible with respect to friction (Reid and 
Bretschneider, 1953). In addition, external wind influences 
are neglected. Kajiura (21964) discusses the partial 
Mes Lection Of water waves passing over a bottom or 
Variable depth. Green's formula provides a very good 
estimate of the transmitted wave amplitude even 17 
partial reflection exists on a sloping bottom (Reid, 
I 5 796 

In the absence of energy losses by friction ard 
۱۳ ۳ 1100 2ه‎ dispersion Dy Cr fraccion and reflection 
the wave energy transmitted between orthogonals (wave 


rays) remains constant. 
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( ንህ 处 | ۷ 
۰۶۴0 ۳۰۴۰٠۷۰ ٣٢ت ۸ر۸‎ p P ٣٣٣٥9 


1. (Mevesnesrserıon 


Gravity wave refraction is analagous to the bendant 


Ene Lays 1n geometrical optics, Since p w 


N 


Speed decreases with decreasing depth (C - (gh)^), tne 


wave rays are refracted towards regions of lesser water 
860+ 


Dns can be put in ت6‎ 62727066706 6166/۶۶ mg 


m) گرم‎ 


wnere ሠ is the angle of the wave ray measured counter- 
clockwise from a fixed reference (the x-axis), s is arc 
Mrstanee on the wave ray and n is distance measurea 
Parallel to the wave crest and to the left of the ray. 
The above relation can be shown to be consistent with 
Fermat's principle (Munk and Arthur, 1952). In rhe spec- 
mp ecco wor Ssrraight and barallelydeptnwuconcoursysrore 


Merece “ne x=axis normal Eo the contours such ina. 


2:280 ones in thts case (l1. I reduces co 
m D 
SORGE, ھ9‎ ኋ sin [J ہے‎ 
dx C dx 
which 1s readily integrated and yields Snell's law- 











pu 


CE) sin. 

Core ጋግ 
elong the wave ray. Note however, that tnis is a mmeceiz! 
case of the more general relation (1.1). 

IN ۱5 (۱۵۵1 و2‎ ርን ፡ጅ.1 l methods arca S nn 
referred to in describing the refraction pattern. The 
two methods most frequently employed are: first, the 
crest method, where successive positions of tne wave 
er esL are drawn by plotting che wave advance Irom point 
me point along the crest; second, the ray method, where 
Pre gor .no00nel as plotted direcely Dy determinirg ai 
snoreward deflection as it crosses successive bottom 
contours. Johnson, O'Brien, @nd Isaacs (1948) provide 
۴۰۰۰۰٠۰٠۶۰۰۶۰۰٠۰ 0ت ت0 667 1ط‎ ۳ 6 7 ፎሮ”  503:5ሀ:(.11::01656፡..)1.፡፡፡፡6፡ ሽብ ጉ፡. 
a detailed large scale chart of the bottom topograpny 
is essential. For the advancages and disadvantages orf 
each method the reader is referred to an excellent Qis- 
bussropD"5yv Dunnham (1950). 

Ne rrer reaction Ol -energy in ۳ 1 
we (Jep o determine by using a gr£ésphical approsen, 
እ. ከ1 ال‎ reason analytical solutions for teje rays neve 
seen qeveloped, where the depth. variation is expressec 
ES uoc 00ب تج‎ 0 (Arthurs 1 7, Srerore prosee 


۰۶۰۰۰۰۰۰۹۰ CISCUSSION of tais par Ul aslana Case 


0-2.0 
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(0 


enesr lorıonsnın between. wave nec M Dd enercv 





reviewed. 


Ir it is assumed trde GREP cotcei Mercy trarnswibctec 


(V 


Derween two orthogonals is constant, convercence of the 
tays denotes a concentration orf energy and a correspond- 
ing increase in wave heights. To z= good approximation, 
Dave height is proportional to che souere root or tre 
energy provided that the waves are not near tne break- 


۳ ¿Tt (Murk and Traylor, 1947). The wave helc c, 


zn 
ጨ 


= 


however, will not usually be the same along a particu 


= 


ED resc due to refraction effects over an 8 ٦ 


bottom. 
The mean wave energy per unit surface area for 


progressive waves is given by 
pe 
a Se 


asco O is the density of sea water, g the acceleretion 


Bue co gravity and H the wave height. The mean weve 


۱2۱ 
O 
pa 
| = 
0 
= 
U} 


PONS EG 1S constant between two wave FCS sc 


ግ ፥ « ۰ 


tnat the relationship for tne wave heignt is given ۷ 


es 
ED 


IN 


58 
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where Ha is the deeo water wave height, Kis che 


peGrractron' factor and X. ج27‎ ٠٣۰ ጎጎ” በጋ 





= 
Ps 


Mme refraction fector 2s ۰)۹ ۶ by 


Mo 4) K = (b E 


= 1/2) 


where b is the distance between adjacent rays measured 
normal to the rays and the subscript l refers to deep 
water 


The shoaling factor can be exoressed bv 


(1.5) K. = (V./V) 


ممه 


where Va is the deep water group speed. In general 
VREA 


where n is given by 


For the wave periods considered in this study n is 


.ታ. ትም l y unity ror all depens and hence oV WX co ILS 
ር እ. 
1.0) K. = (h/h)? æ 


G ان‎ ncident wave O unit اد ها تا‎ oO = 
تحت‎ As mand) (156) thatthe relai o I GEC EES 


emo@itude is given by 


(1.7) je Kn E š 





is 
2-٥7 0۲016 SG C iphe 


Fermat's principle o- Base come 25 ፡መሰለ لا‎ nn 
ec IGA Tt can be steteo thus: (Ene ropa s NN . 
light always takes place in such a way that the actual 
optical path (e.g. length of geometric path multiplied 
by index of refraction of the medium) is an extreme 
value compared with all other oaths whicn do not follow 
meee few ©f Optics, In @ weve refraction analogy, we allow 
pne soHeOreys to correspond to ortroconals AF a gravii 
wave refraction diagram. 

The application of Fermat's least time concesc to 
N TION وت لد 1ے ۲۱۹۷66 طط‎ e rC olar ۴0 0 =. 
discussed by Arthur (1946). The travel time along a path 


一 - 


Bon a point A tO a DOINE 5 Of a refracted ray par .. 
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Dy 


mud) 13 = Gere 9 
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1 
A - ۹ ۰ ' 
where C s (gh)^ along the peth and ds is the arc length 


ry 


long the path. Consider depth h = hlrJ] only; chef fos 
I to be a minimum recuires an r(€) such chat $I = Ue 
Bor EE, the depth is assumed constant and tne or- 


Mreegome SG are streignt lines.) For © = ሚሚ 


= the Jm NES 
O 2 in 


are curved. Using the coordinates es shown in 


۱ E 
Le 


By (221) takes the form 





56 
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The refraction factor is given by the equation 
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A calculation summary for the circular island 
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The phase lag on the island shore increases with 
سے‎ (see Figure 11). The range of phase lag at tne 
E rand shore varies from about l radian for waves orf 
12 minute period to about 6.5 radians for waves ol 2 
Minute periodes TInus the phase lag of waves in tne lee 
# tne island varies considerably with period, while 


the amplitude is unaffected by the period (within the 


tsunami range). 
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